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The  research  conducted  under  Grant  Number  AFOSR- 83-0210  during  the  period  15 
July  1983  to  14  July  1984  is  reported  on  in  G>]»  [5]  and  [6]. 

The  goal  of  this  project  is  to  develop  methods  to  solve  efficiently  constrained 
minimization  problems  that  have  problem  functions  that  are  not  everywhere  differentiable. 
Such  difficult  problems  occur  in  practice  when  decomposition,  nested  dissection, 
relaxation,  duality  and/or  exact  penalty  techniques  are  applied  to  large  or  complicated 
nonlinear  programming  problems  in  order  to  convert  them  to  a  sequence  of  smaller  or 
less  complex  problems.  Being  able  to  use  these  techniques  that  often  give  implicitly 
defined  problem  functions  gives  a  user  flexibility  in  modeling  a  problem  for  solution 
and  the  ability  to  exploit  parallel  processing  in  computation. 

In  (5 ]  the  principal  investigator  developed  a  reliable  and  efficient  method  for 
solving  single-variable  optimization  problems  that  may  have  constraints  and  may  have 
problem  functions  that  are  nonsmooth  and/or  nonconvex.  The  method  converges  to  stationary 
points  of  problems  on  which  secant  and/or  cutting  plane  methods  may  fail  and  it 
converges  rapidly  on  some  problems  on  which  these  other  methods  do  not  exhibit  super- 
linear  convergence.  The  new  safeguarding  development  that  insures  convergence  via  a 
sequence  of  intervals  with  decreasing  lengths  does  not  detract  from  better  than  linear 
convergence  and  fits  in  well  with  numerical  practice,  because  it  keeps  apart  points 
whose  difference  appears  in  the  denominator  of  a  divided  difference  formula  for 
second  derivative  approximation.  Our  new  automatic  penalty  technique  for  handling  a 
constraint  Is  also  numerically  sound,  because  It  is  independent  of  constraint  scaling 
and  numerically  well -conditioned  under  a  mild  constraint  qualification. 

This  new  single  variable  method  can  be  used  to  solve  line  search  problems  in 
multivariable  optimization  where  a  line  search  function  value  decrease  requirement  can 
be  Imposed  via  a  constraint.  Also,  It  can  be  used  to  solve  a  root  (or  zero)  finding 
problem  on  an  Interval  by  minimizing  the  absolute  value  of  the  appropriate  function 
subject  to  a  constraint  defining  the  Interval. 

To  show  the  efficacy  of  our  single  variable  method  and  corresponding  FORTRAN^de 
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we  solved  a  practical  single  resource  allocation  problem  [1]  with  five  bounded  decision 
variables  via  a  dual  or  min-max  technique.  This  application  used  our  subroutine  in  a 
nested  manner,  i.e.  a  single  variable  "outer"  problem  was  solved  where  each  function 
evaluation  involved  solving  a  five  variable  "inner"  Lagrangian  problem  that  separated 
into  five  independent  single  variable  problems.  This  nested  application  of  our  algorithm 
used  somethat  less  function  evaluations  that  the  number  used  by  Powell's  sequential 
quadratic  programming  code  VMCWD  [7]  applied  to  the  primal  problem.  If  one  has  the 
facility  for  parallel  processing  then  our  approach  can  be  made  even  more  efficient  by 
solving  the  five  independent  subproblems  in  parallel.  We  reported  these  results  at 
the  SIAM  Conference  on  Numerical  Optimization  held  at  Boulder,  CO  in  June  1984.  We 
have  been  invited  to  write  a  code  description  and  user's  guide  for  our  FORTRAN  imple¬ 
mentation  for  presentation  at  the  NATO  Advanced  Study  Institute  on  Computational  Mathe¬ 
matical  Programming  to  be  held  at  Bad  Windshelm,  West  Germany  in  July  1984.  After 
completing  this  we  intend  to  revise  and  Improve  our  earlier  paper  on  a  BASIC  imple¬ 
mentation  [2]  so  that  challenging  problems  can  be  solved  on  micro-  as  well  as  on  more 
sophisticated  computers. 

The  principal  investigator  and  his  Ph.D.  student  Nilama  Gupta  are  making  progress 
on  generalizing  the  above  research  to  the  multivariable  case.  Part  of  this  work  for 
the  unconstrained  case  was  presented  at  the  11th  IFIP  Conference  on  System  Modelling 
and  Optimization  held  at  Copenhagen,  Denmark  in  July  1983,  and  will  appear  in  [6]. 

This  research  Involves  developing  a  safeguarded  piecewise  quadratic  approximation 
technique  that  approximates  higher  than  first  order  Information  and  leads  to  a  new 
concept  of  better  than  linear  convergence.  This  method  is  being  designed  to  converge 
to  stationary  points  of  problems  with  semi  smooth  functions  (a  concept  developed  under 
Grant  Number  AF0SR-74-2695  [3]).  It  should  converge  rapidly  under  stronger  assumptions 
including  an  underlying  piecewise  C^  structure. 

The  proposed  method  requires  the  solution  of  a  certain  Indefinite  quadratic  pro¬ 
gramming  subproblem,  so  our  Ph.D.  student  Amal  Al-Saket  is  developing  a  numerically 


reliable  method  to  solve  such  problems  (and  generalizations 


thereof)  based  upon  our 


earlier  successful  efforts  on  the  constrained  least  squares 


problem  [4], 
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